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The light-front dynamics is an efficient approach to study of field theory and of rela- 
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■ tivistic composite systems (nuclei at relativistic relative nucleon momenta, hadrons in the 

quark models). The explicitly covariant version of this approach is briefly reviewed and 
\ illustrated by some applications. 

1 Introduction 

o 

The recent experimental data on deuteron structure and heavier nuclei concern the relativistic 
domain of the nucleon momenta. Their successful theoretical description was achieved |]. Q 
in the framework of light-front dynamics (LFD). This approach deals with the wave function 
$_i ■ denned on the four-dimensional space-time plane given by equation t + z = 0. 



In general, the relativistic state vector is defined on a space-like surface. The plane 
t = const, associated with the non- relativistic quantum mechanics, is a particular case of this 
surface. The Schrodinger equation describes the dynamical time-evolution of the wave function 
ijj(r,t), or, in other words, from one space- like plane t = const to other one. This dynamics 
was called by Dirac the "instant form" ||. Light-front dynamics, also called "front form", is 
associated with the plane t + z = 0. The latter is considered as a limiting case of a space-like 
plane. For completeness we mention the third form proposed by Dirac: the "point form", 
associated with the hyperboloid x 2 = t 2 — r 2 = const. 

This paper is a brief review of the explicitly covariant version of light-front dynamics. We 
explain its relation to the Bethe-Salpeter approach and illustrate it by few examples. Namely, 
we give the wave function of two relativistic scalar particles interacting by scalar exchange, dis- 
cuss the relativistic deuteron wave function and its electromagnetic form factors in comparison 
to CEBAF/TJNAF data and present our new results concerning the system of three relativistic 
bosons with zero-range interaction. 

*e-mail: karmanov@sci.lebedev.ru, karmanov@isn.in2p3.fr 
^e-mail: carbonel@isn.in2p3.fr 
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1.1 Why light-front dynamics? 

The theoretically self consistent description of the relativistic composite system is based on 
the quantum field theory. The state vector of a composite system and even of the vacuum 
has a very complicated structure. It is usually described as a superposition of bare quanta, 
corresponding to non- interacting fields. If we "switch off" the interaction between fields, the 
number of particles is conserved. As soon as we take into account the interaction, the state 
vector, is in general a superposition of states with different numbers of particles. 

When the interaction is weak, like in the case of quantum electrodynamics, it does not 
change the state vector too much and the "dressed" electron differs from the bare one mainly 
by small admixture of photon. 

The situation is drastically different, when the interaction is strong. For example, the 
proton is usually treated as being a three quark system, but these quarks are not the same than 
those appearing in the initial Lagrangian of Quantum Chromodynamics. They are constituent 
quarks, which, in their turn, consist of bare quarks and gluons. The proton state vector is a 
huge superposition of bare fields and has not yet been calculated from the QCD first principles. 

One should emphasize that not only the proton state, but also the state without physical 
particles - the vacuum state - is a complicated superposition of the bare particles, or, in other 
words, of vacuum fluctuations. This description of emptiness in terms of the very complicated 
conglomerates of particles seems unnatural. It would be much better to use an approach, in 
which the vacuum is indeed nothing but emptiness. Simplifying the vacuum wave function, 
we simplify the wave function of the dressed states as well, eliminating from them, like in the 
vacuum wave function, the field fluctuations. After doing that one can study the real, physical 
structure of particles and nuclei in a more higienic way. 

In the light-front dynamics the vacuum is nothing but emptiness and here is the main 
advantage of this approach. 

Qualitatively it can be understood, by noticing that the front form corresponds to the 
instant form, studied from a reference system moving with velocity close to the speed of light. 
Indeed, an observer moving with a velocity v along z-axis describes a physical process in his 
coordinates (t 1 , x', y', z'), which are related to the laboratory ones by the Lorentz transformation: 

. z — vt . t + zv/c 2 . . 
z = i t = ^=, x = x, y =y. 

yjl - V 2 /C 2 sjl - V 2 /C 2 

The plane t' — in the moving system corresponds to t + zv/c 2 = in the laboratory 
coordinates. In the limiting case, when v — > c, we get the light-front plane t' oc z + = t + z/c = 0. 
So, the front form corresponds to the normal form if — 0, but studied from the system of 
reference moving with the limiting speed c called since Weinberg Q the "infinite momentum 
frame" . 

Consider now the fluctuation creating three particles from vacuum (e.g., in the y9 3 -theory) . 
The fluctuation with energy AE = +£% 2 Jr£ k 3 ma y occur during the time At ~ h/AE (here 

ej: = \Jk 2 -\- m 2 ). In the infinite momentum frame the momenta fcj and energies eg. of any 
particle increase and AE tends to infinity. Therefore, the fluctuation time At tends to zero. 
The contribution of this fluctuation to the vacuum wave function disappears. Due to the time 
dilation, all physical processes are delayed, and fluctuations have no time to occur. This means 
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that in a thought experiment in the infinite momentum frame we study the particles prepared 
"far in advance", not spoiled by the vacuum fluctuations. From practical point of view, it is 
enough to work in a usual reference frame, but with the state vector defined on the light-front 
plane. 

A disadvantage of this approach is the distinguished role of £-axis relative to the x, y-axes. 
Because of that, the theory is not explicitly covariant. To overcome it, in the paper || (see 
for review || [7|) the explicitly covariant version of light-front dynamics (CLFD) was proposed, 
which deals with the state vector defined on the light-front plane u-x = 0, where uj = (uj ,uj), 
uj 2 = uj 2 — uj 2 = 0. The theory becomes explicitly covariant, whereas the dependence of the 
wave function on the light-front plane is parametrized by means of the four-vector uj. This is 
very convenient from the point of view of applications. In particular the standard approach is 
recovered when to = (1, 0, 0, —1). 



2 Wave function 

The wave functions studied in LFD are the Fock components of the state vector, that is, the 
coefficients in an expansion of the state vector \p > with respect to the basis of free fields: 

\p) = (2tt) 3 / 2 U 2 {k l M^r)a){k 1 )a\k 2 )\Q) 



x 5^( kl + k 2 -p- ^)2(^)dr w ^ 2 _ (27r) 3/?^ + • • • ■ (1) 

The dots include the higher Fock states. For simplicity, we omit the spin indices. By con- 
struction, all the four-momenta are on the corresponding mass shells: k 2 = m 2 , p 2 = M 2 , 
(lut) 2 = 0. 

We emphasize the presence in (|IJ) of the delta-function. It provides the conservation law: 

h + k 2 = p + ujt, (2) 

where r is a scalar variable responsible for the off-energy shell continuation. The four-momentum 
ujt is associated, for convenience, with a Active particle called " spurion" . We would like to em- 
phasize at this point that the use of " spurion" does not imply the use of any unphysical degres 
of freedom in the theory but only a confortable way of parametirsing the off-shell effects. As 
mentioned above, the standard LFD corresponds to uj = (1, 0, 0, —1), or, in the light-front coor- 
dinates, uj + = ujq+u z — 0, cl>j_ ee u x>y = 0, uj- = luq—uj z = 2. In this case, equation (fj) gives the 
standard conservation laws for the (_L, +)-components of the momenta (which does not include 
any r), but does not constrain the minus-components. Their difference k\_ + k 2 _ — p_ = 2r 
determines the value of r. Because of non-conservation of the minus- components, the wave 
function is "off-energy shell". 

The term "off-energy shell" is borrowed from the old fashioned perturbation theory, where 
it means that for an amplitude which is an internal part of a bigger diagram, there is no 
conservation law for the energies of the incoming and outgoing particles. Note that any bound 
system is always off-energy shell, since the sum of the constitutent energies le^ (in the c.m.- 
frame of two-body system) is not equal to mass M of the bound system but larger: 2e^ > 2m > 
M. We emphasize the difference between the off-energy shell and off-mass shell amplitudes: 
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in the latters the four-momenta squared are not equal to the masses squared, but satisfy the 
conservation law for all the four components. An example of off-mass shell amplitude is the 
Bethe-Salpeter function (see below sect. ^). 

As mentioned, the wave function, defined on the light-front plane, depends on orientation 
of this plane. From ([I]) we see that this dependence is given by the four- vector uj. It is an 
important property of any Fock component. At the same time, the observables, like form 
factors, don't depend on u. In principle, the form facors, being scalars, could depend on the 
scalar products F = F(q 2 , uj-p, uj-p'), where q 2 = (p — p') 2 . Since u determines a direction only, 
the theory is invariant relative to the replacement uj — > au, where a is a number. Therefore, 
the extra scalar products should enter in the form of ratio: uo-p/uj-p' . To factorize the e.m.- 
amplitude, one should impose the condition uj-q = (equivalent to q + = in the standard 
approach, see J7| for details). With this condition we get uj-p/uj-p' = 1, so the form factor does 
not depend on uj. 

For the states with spins 1/2 and 1, decomposing the e.m. vertex in the form factors, one 
can construct not only scalar products and not only usual spin structures (two for spin 1/2 and 
three for spin 1), but also tu-dependent spin structures. The scalar coefficients at the front of 
them are extra, non-physical form factors. One should separate the physical form factors from 
non-physical ones. Corresponding formulas extracting the physical form factors from full e.m. 
vertex were found in || [!|. 

— # 

Parametrized in terms of the relative momentum k, the two-body wave function reads: 

ip = ip(k,n), (3) 

where k = k\ in the system where k\ + ki = and n is the direction of uj in this system. We 
see that the relativistic light-front wave function depends not only on the relative momentum 
k but on another variable - the unit vector n. 

It is also useful to introduce another set of variables 

x = uj-ki/uj-p , R = ki — xp (4) 

and represent the spatial part of the four- vector R as R — R» + R±, where R\\ is parallel to u 
and -R_|_ is orthogonal to uj. Since Ruj = R^ujq — Rpjj = by definition of R, it follows that 
Ro = \R\\\, and, hence, R\ = —R 2 is invariant. Therefore, R\ and x can be chosen as two 
scalar arguments of the wave function: 

i> = i>(R 2 ± ,x) . (5) 
Using the definitions of the variables R\ and x, we can readily relate them to k 2 and n-k: 

R 2 ± = k 2 -(n-k) 2 , ar = I M - 

The variables introduced above can be easily generalized to the case of different masses and an 
arbitrary number of particles |§. 




4 



3 Equation 



The equation for the wave function is shown graphically in fig. |l| The dashed line corresponds 
to "spurion". In the case of a scalar particle, in terms of the variables k,n, this equation has 
the following form: 

(A(P + m 2 ) - M 2 ) n) = -^3 f i){k\ n)V(k', k, n, M 2 )— . (7) 
The kernel of eq. (0) depends on the vector variable n and on the bound state mass M 2 . This 




Figure V. Equation for the two-body wave function. 

dependence is associated with the retardation of the interaction. 

In the non-relativistic limit the light front wave function turns into the ordinary equal- 
time wave function. This limit can be formally obtained at c — ► oo. Then the light-front plane 
t + z/c = turns into the plane t — 0. Therefore the equation (0) turns into the Schrodinger 
equation in momentum space, the kernel V being the Fourier-transform of non-relativistic 
potential and the wave function no longer depends on n. 



4 Simple example 

For a given Lagrangian, the amplitudes are calculated by rules of special graph technique, 



developed in and adjusted to the light-front dynamics in ||. The diagrams correspond 
to the time-ordered graphs. The amplitudes on energy-shell coincide with usual Feynman 
amplitudes on-mass shell. Off the energy and mass shells, these amplitudes differ from each 
other. 

Consider exchange in the t channel by a scalar particle of mass /i between two scalar 
particles. The amplitude is given by two time-ordered diagrams shown in fig. |2]. They determine, 
in the ladder approximation, the kernel of the equation (|7|). The external spurion lines indicate 
that the amplitude is off-energy shell. 

According to the light-front graph technique, the amplitude has the form: 

K = g 2 6 (ujjh - k[)) + g 2 6 (u{k[ - h)) 



- (h - k[) 2 + 2tw(A; 1 - k[) -te /i 2 - (k[ - h) 2 + 2r'uj{k' 1 - h) - ie 
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Figure 2: Exchange by a particle in t-channel. 



The two items in (||) correspond to the two diagrams of fig. |2|. They cannot be non-zero 
simultaneously. The off-energy shell terms, depending on u>, appear automatically and they 
are explicitly given in (||). On the energy shell, i.e. for both r = t' = 0, the expression for the 
kernel becomes identical to the Feynman amplitude: 



o) 



IE 



(9) 



In the case of /i = we obtain the Wick-Cutkosky model pT| , but treated in the framework 
of CLFD |Tj|. Going over from the kernel K to V — — /C/(4m 2 ), introducing the constant 
a = g 2 /(16iim 2 ), and expressing (Q) by means of the initial and final relative momenta k, k', 
we obtain 1X211 : 



V 



4-rra 



ik' -kf -{n-k'){n-k) 



+ et + ei 



2 ; 



n-k ' 



£k> 



n-k 



(10) 



The dependence on the light-front orientation manifests itself in ( |i0|) through the vector pa- 
rameter n. 

The solution of the equation (|7|) with the kernel (JlQ|), in the limit of small binding energy, 
has the form: 

8^nmn 5 / 2 



i/j(k, n) 



{k 2 + K 2 ) 2 



1 + 



\n-k\ 



11) 



In non-relativistic limit k, k! <C m the kernel (|l0|) turn into the Coulombien one: V ~ 
— 4ira/(k ' — k ) 2 (that in the coordinate space corresponds to V(r) = —a/r). The n-dependent 
factor in the wave function flTTp disappears and the latter turns into the Coulombien wave 
function too. 

We emphasize that the dependence of the wave function ( |TTD on n does not mean any 
violation of the rotational invariance, since this function is a scalar, as should be for the zero 
angular momentum. However, for non-zero angular momentum, because of Fock space trun- 
cation, one can find dependence of energy on the projection of the angular momentum on the 
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n-direction. This fact indeed means a violation of rotational invariance. A solution of this 



problem is given in 13 



The solutions of equation (|7|) with the kernel ( JT0| ) for < fi < m were studied in [|14]]. It 
was found that for the exchange mass [i comparable with the constituent mass m the relativistic 
effects in the binding energy are very important even for vanishingly small ("nonrelativistic") 
binding energy 



5 Relation with the Bethe-Salpeter function 

It is instructive to compare the solution ([TT| ) with one found using the Bethe-Salpeter function 



p"5fl . The latter is defined as (see for review |T6|] ): 

$(x u x 2l p) = (OlTC^rO^))^). (12) 

Here ip(x) denotes the field operator in Heisenberg representation. The function fli~2| ) depends 
on two four- vectors xi,x 2 , they include two times ti,t 2 - Though the Bethe-Salpeter function 
satisfies the normalization condition, allowing to find the normalization coefficient, it has no 
any probabilistic interpretation. 

In the momentum space the Bethe-Salpeter function looks as: $ = &(qi,q 2 ,p) and is 
defined through the Fourier transform: 

$(xi,x 2 ,p) = (27r)~ 3/2 exp [—ip-(xx + x 2 )/2] <f>(x,p) , x = x x — x 2 , 

Q2,p) = <&(q,p) = J &(x,p) exp(iq-x)d A x , (13) 

where q = (q± — q 2 )/2, p = qi + q2, qi and q 2 are off-mass shell four-vectors: q\ ^ m 2 , q\ ^ m 2 . 

In LFD the Heisenberg operators coincide with the Schrodinger ones on the light-front 
plane (like in the ordinary formulation of field theory the Heisenberg and Schrodinger operators 
are identical for t = 0). Therefore, projecting the Bethe-Salpeter function on the light-front, 
we can relate it to the light-front wave function. In the momentum space this relation reads 



1>(k 1 ,k 2 , P ,UT)= t^ 1 )^ 2 ) / +00 $( gi = k 1 -u}T/2 + uP,q 2 = k 2 -UT/2-uj/3,p)d/3 (14) 

ir{U!-p) J-oo 

where &(qi,q 2 ) is the Bethe-Salpeter function fll3|). The argument p in ([14]) is related to the 
on-shell momenta ki,k 2 as p = k\ + k 2 — out, that follows from the off- mass shell relation 
P = qi + <?2- 

The equation (0) makes the link between the Bethe-Salpeter function $ and the light- 
front wave function ip. It should be noticed however that the wave function (|14D is not neces- 
sarily an exact solution of eq. (0), since, as a rule, different approximations are made for the 
Bethe-Salpeter kernel and for the light-front one. In the ladder approximation, for example, 
the second iteration of the Bethe-Salpeter amplitude contains the box diagram, including the 
time-ordered diagram with two exchanged particles in the intermediate state. This contribution 
is absent in the light-front ladder kernel and in its higer order iterations. 
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We can check now that in the Wick-Cutkosky model the relation ( 111 ) reproduces the wave 
function (|TTD . In this model, the exact expression for the Bethe-Salpeter function has the form 



$(g,p) = — ic 



m 2 - ^M 2 - q 2 



m 2 — (-p + q) 2 — iO ) ( m 2 — {—p — q) 2 — /() 



-i 



(15) 



where c = 2 5 V ttttik 5 with k = */m|e&| = ma/2. Substituting ( |I5"D in (0), we indeed reproduce 



the expression (|TT|) for the light-front wave function. 

The e.m. form factors calculated through the light-front wave function and through the 
Bethe-Salpeter one were compared in M. They coincide with each other with high accuracy. 
This means that when both approaches take into account the same relativistic dynamics, they 
give close results. 

However, an important difference between them is in the following. Bethe-Salpeter func- 
tion (|12|) is a formally defined object. Whereas, the light-front wave function, being the co- 
efficient in the decomposition ([]]), according to general principles of quantum theory, is the 
probability amplitude. It is the same object as nonrelativistic wave function and it has the 
same physical meaning. Comparing light-front wave function of a relativistic system with its 
non-relativistic counterpart, we describe the relativistic effects in a language of physical phe- 
nomena. This important point is an advantage of LFD. Smooth matching of non-relativistic 
and relativistic wave functions is also very convenient in practical calculations. 



6 Deuteron wave function and form factors 

The deuteron wave function is determined by six invariant functions, depending on two scalar 
variables. This number is the dimension of the matrix depending on the spin projections 
of the deuteron and two nucleons, divided by the factor 2 due to the parity conservation: 
N = 3x2x2/2 = 6. This wave function reads: 



^(M) = ^^^n^" 3 2 ( ( )_ 1 



+ h— (3k(n-a) + 3n(k-a) - 2{hn)a) + } ; 



k x n 



(16) 



+ [[k xn}xo 



In the relativistic one boson exchange model it was calculated in |17[|. The results are displayed 



in fig. One can see that the function f$, of relativistic origin, is very important: it dominates 



at k > 500 MeV/c. In [|18l it was shown that this component dominates since it automatically 
takes into account the isovector exchange currents (so called pair terms). They should not be 
added separately. In nonrelativistic limit the functions /3_6 become negligible and only two 
first functions /i_2, corresponding to usual S- and D-waves, survive. 

This wave function was used in fl], to calculate the deuteron electromagnetic form 
factors. No any parameters were fitted. It turned out that the calculated polarization observable 
£20, fig- H, as well as the charge and quadrupole form factors, fig. ||, coincide with rather precise 
experimental data obtained recently at CEBAF/TJNAF P"5[ . 
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Figure 3: Deuteron relativistic components f\ — / 6 as a function of k, at nk/k = 0.6. 



7 Three-boson bound states with zero-range interaction 

The most recent application of CLFD is the relativistic system of three-bosons with zero range 
interaction. We will describe it in more detail. 

Zero-range two-body interaction provides an important limiting case which qualitatively 
reflects characteristic properties of nuclear and atomic few-body systems. In the nonrelativistic 
three-body system it generates the Thomas collapse J21]. The latter means that the three-body 
binding energy tends to — oo, when the interaction radius tends to zero. 

When the binding energy or the exchanged particle mass is not negligible in comparison 
to the constituent masses, the nonrelativistic treatment becomes invalid and must be replaced 
by a relativistic one. Two-body calculations show that in the scalar case, relativistic effects are 
repulsive (see e.g. [|14||). Relativistic three-body calculations with zero-range interaction have 
been performed in a minimal relativistic model |22| and in the framework of the Light-Front 
Dynamics p3 |. In these works it was concluded that, due to relativistic repulsion, the three- 
body binding energy remains finite and the Thomas collapse is consequently avoided. However, 
in the papers |2^, ^3] a cutoff has been introduced implicitly. Being interested in studing the 
zero-range interaction as an important limiting case, we eliminate this cutoff. 

The three-body equation is represented graphically in figure || It concerns the vertex 
function T, related to the wave function ip in the standard way: 

ip[kx,h,h,p,ujT) = M2 _ — 2 , M ={k 1 + k 2 + k 3 ) ={jp + ut). 

Like for the two-body system, all four-momenta are on the corresponding mass shells: k\ = 
m 2 ,p 2 = Mf , (wr) 2 = and satisfy the conservation law k\ + k 2 + k 3 = p + lot involving lot. 
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Figure 4: Polarization observable £20 i n &d scattering, taken from |20| 
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Figure 5: Charge and quadrupole deuteron form factors, taken from | 
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Figure 6: Three-body equation for the vertex function T. 



Applying to figure ^| the covariant light-front graph techniques [0] , we find: 

A 



k 2 , k 3 ,p, lot) 



(27T) 



: . , r(k[,k' 2 ,k 3 ,p,ujT') 



fin- 1 rl 3 h' rl 3 h' 

x5^\k[ + k' 2 -Lor' -h-k^ + ior)^-^^ + (23)1 + (31)2, 



t' 2e k > 2s k 



(17) 



where Sk = sjrn? + k 2 . For zero-range forces, the interaction kernel is replaced in fllTD by a 
constant A. The contribution of interacting pair 12 is explicitly written while the contributions 
of the remaining pairs are denoted by (23)1 + (31)2. 

The Faddeev amplitudes are introduced in the standard way: 

r(i, 2, 3) = r 12 (i, 2, 3) + r 23 (i, 2, 3) + r 3 i(i, 2, 3), 

and equation (|I7]) is equivalent to a system of three coupled equations. With the symmetry 
relations r 2 3(l,2,3) = r 12 (2,3, 1) and r 3 i(l,2,3) = Ti 2 (3, 1,2), the system is reduced to a 
single equation for one of the amplitudes, say T 12 . 

We use the tree-body version of the variables defined in eq. In general, r i2 depends 
on all variables (Rij_, Xj), constrained by the relations Ru_ + R 2 j_ + R 3 ± = 0, X\ + x 2 + x 3 = 1, 
but for contact kernel it depends only on (-R3_l,x 3 ) p3| . 

The value of A is found analytically by solving the two-body problem with the same zero- 
range interaction under the condition that the two-body bound state mass has a fixed value M 2 . 
In this way, the constant A disappears from the kernel, but the two-body scattering amplitudes 
appears, as in any Faddeev equations. For the zero-range interaction, this amplitude does not 
depend on the relative momentum, but depends on the off-energy shell variables. Therefore it 
is extracted from the integral. Finally, the equation for the Faddeev amplitude reads: 

r 12 (R' ± , x'(l - x)) d 2 R' ± 



T 12 (Rj 



,x 



F(M 



12; 



dx' 



(2nyJo Jo (R' ± - X 'R ± ) 2 + m i- x '(i- X ')M : 
where F(Mi 2 ) is the two-body off-shell scattering amplitude: 

arctany^ 



F(M_ 



12 



2y 



Ah 



1 1 + VMi2 

log _ r ia 

1 - VM 12 



Vm 2 



SIT' 



arctan Um 12 arctan t/M 2 



VMy. 



Vm 2 



M\ 2 < 0, 



< M\ 2 < Am 2 , 



VMv 



VMx 



Mi 
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and similarly for y M2 = —j=L==. The variable M 2 2 = (k[ + k 2 — cor') 2 = (p — k 3 ) 2 is the 
off-shell mass of the two-body subsystem. It is expressed through Mf , R\, x as 

M ?2 ^(l-,)M|- fl ' + ' 1 -^ . (19) 

The three-body mass M 3 enters in the equation ([T^) through variable M 2 2 . 
By replacing x'{l — x) — > x', equation (|l"8D can be transformed to 

r 12 (^,x) = F(M 12 )-— 3/ ^ / _|i^_r 12 (^,x'), (20) 



(27r) 3 Jo ar'(l - x - x') Jo .M -Mi 



where 



^ /2 _ R'l + m 2 + R\ + m 2 + (i2' ± + # ± ) 2 + m 2 



This equation is the same than equation (11) from ]23| except for the integration limits 



of (R'±,x') variables. In []23] the integration limits follow from the condition M\ 2 > 0. They 
read 

[. . .} dx' / X [. . .] d 2 R' ± (21) 



with k™ ax = y (1 — x')(M 2 x' — m 2 ) and introduce a lower bound on the three-body mass 

M 3 > y/2m. The same condition, though in a different relativistic approach, was used in |[22j| . 
The integration limits in (|2l| ) restrict the arguments of T 12 to the domain 

2 9 

77? 777 

JWf- " M 2 ' 

and can be considered as a method of regularization. In this case, one no longer deals with the 
zero-range forces. 

The results of solving equation ( |18| ) are presented in what follows. Calculations were 
carried out with constituent mass m = 1 and correspond to the ground state. We represent 
in fig. |^a the three-body bound state mass M 3 as a function of the two-body one M 2 (solid 
line) together with the dissociation limit M 3 = M 2 + m. The zero two-body binding limit 
B2 = 2m — M 2 — ^ is magnified in fig. |7|b. In this limit the three-boson system has a binding 
energy B s m 0.012. 

When M 2 decreases, the three-body mass M3 decreases very quickly and vanishes at 
the two-body mass value M 2 = M c pa 1.43. Whereas the meaning of collapse as used in the 
Thomas paper implies unbounded nonrelativistic binding energies and cannot be used here, the 
zero bound state mass M 3 = constitutes its relativistic counterpart. Indeed, for two-body 
masses below the critical value M c , the three-body system no longer exists. 

The results corresponding to integration limits ( 21 ) are included in fig. [7]a (dash line) for 



comparison. Values given in [fj3] were not fully converged |25|. They have been corrected in 
23J and are indicated by dots. In both cases the repulsive relativistic effects produce a natural 



cutoff in equation (|IBD , leading to a finite spectrum and - in the Thomas sense - an absence of 



collapse, like it was already found in f22j| . However, except in the zero binding limit, solid and 



dash curves strongly differ from each other. 
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Figure 7: (a) Three-body bound state mass M 3 versus the two-body one M 2 (solid line). Results 
obtained with integration limits ([H]) are in dash line. Dots values are taken from |24[. (b) Zoom 
of the zero two-body binding region (M 2 — > 2m, _B 2 — > 0) displaying solid line only. 



M, 



-10 - 



Figure 8: Three-body bound state mass squared M| versus M 2 . 



We would like to remark that for M 2 < M c , equation (|18D posses square integrable 



solutions with negative values of M|. They have no physical meaning but Mf remains finite 
in all the two-body mass range M 2 G [0,2]. The results of Mf are given in figure [| When 
M 2 -> 0, M| tends to w -11.6. 
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From the above results we conclude that the three-body bound state exists for two-body 
mass values in the range M c pa 1.43m < M 2 < 2m. At the zero two-body binding limit, the 
three-body binding energy is B3 pa 0.012 m. The Thomas collapse is avoided in the sense that 



three-body mass M3 is finite, in agreement with [53, 23|. However, another kind of catastrophe 



happens. Removing infinite binding energies, the relativistic dynamics generates zero three- 
body mass M 3 at a critical value M 2 = M c . For stronger interaction, i.e. when < M 2 < M c , 
there are no physical solutions with real value of M 3 . In this domain, Mf becomes negative and 
the three-boson system cannot be described by zero range forces, as it happens in nonrelativistic 
dynamics. This fact can be interpreted as a relativistic collapse. 



8 Conclusion 



We have seen that the explicitly covariant light-front dynamics provides a very satisfactory 
framework, both from a theoretical and practical points of view, for describing relativistic 
nucleons in deuteron and other nuclei. The predicted deuteron electromagnetic form factors, 
calculated without any fitting parameters, are in good agreement with experimental data. 

This approach is also successfully when applied to many other problems and a consequent 
effort covering several directions of research is actually being performed. Here is a list of current 
applications. 

• Deuteron wave function [17] and form factors [[!]. [7|. 

• Nucleon momentum distributions in nuclei ||. 



Leptonic decay of heavy quarkonia |26| . 
Fock-space expansion p?| . 

Nucleon wave function PS[ and neutron charge radius in quark model. 



Two-fermion relativistic bound states |3l], |3TJ . 
Perturbative [32] and non perturbative [Q renormalization. 
Three-boson relativistic bound states p4 [. 



We conclude that the explicitly covariant light-front dynamics is a rather efficient approach to 
relativistic few-body systems and to field theory. 

This work is supported by the French- Russian PICS and RFBR grants Nos. 1172 and 
01-02-22002. 
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